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We investigate the enhanced microwave transmission through the array of metallic coaxial annular 
apertures (MCAAs) experimentally and theoretically. The even-mode and the odd-mode surface 
resonances are clarified from the spatial field distributions and the dispersion diagram. The impact of 
local resonance is thoroughly embodied in the even-mode surface resonant states, while the odd-mode 
■ surface resonances are scaled by periodicity, invariant to different local geometry of the unit cell, and 

invisible in measurements. The enhanced transmission is the collective selections on the interplay 
between the local resonances and the evanescent Bloch wave channels on the surface. Transmission 
. measurements for different inner diameter of the apertures show that the transmissivity extrema 

. ' with respect to the specific angles precisely correspond to the degenerate points in the dispersion 

p^*, diagram of surface resonances. 



^ '■ I. INTRODUCTION 

Extraordinary optical transmission (EOT) through a periodical array of subwavelength nano-holes in metal fLlms[l[ 
• t-H , provides a new way to manipulate the light waves in subwavelength division. The origin of EOT was attributed to the 
O ,■ resonant tunneling of surface plasmon polaritons (SPPs) 2\. Later studies found that the enhanced transmission can 
also be achieved in TerahertzQ or microwave regionQ where the metal surface does not support SPPs but surface 
^2 ■ resonances @. The high transmission peak through more complicated apertures, such as coaxial apertures [H, B, 
rectangular holes[l,Q, etc., tends to be red-shifted or even scaled by waveguide resonance instead of the periodicity [l^. 
Polarization-dependent transmission spectra as well as the optical activities are observed when the array of large-sized 
£>V elliptal holes is adopted[ll|. These studies strongly indicate that the local resonance is crucial and dominant to the 
i-C , formation of EOT under specific circumstances. 

In this paper, we adopt the modal expansion method (MEM) to calculate the enhanced microwave transmission 
through the array of metallic coaxial annular apertures (MCAAs) as well as the dispersion of surface resonances. In 
y—i ■ microwave region, metals are perfect electric conductors (PECs), which enable us to easily characterize the spatial dis- 
persion from the periodic structure. The even-mode and the odd-mode surface resonances arc clarified by calculations 
on the transmission spectra, the spatial field distributions, and the dispersion diagram. The impact of local resonance 
is thoroughly embodied in the even-mode surface resonant states, while the odd-mode surface resonances are scaled 
by periodicity, invariant to different local geometry of the unit cell, and invisible in measurements. And the enhanced 
transmission comes from collective selections on the interplay between the local resonances and the evanescent Bloch 
wave channels on the surface. Transmissivity extrema with respect to the specific incident angles precisely correspond 
to the degenerate points in the dispersion diagram of surface resonances. The results are heuristic to the enhanced 
transmission in the optical region as well. 
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II. SAMPLE DESCRIPTION AND EXPERIMENTAL SETUP 



& : 

Our samples were fabricated on lmxlm copper slabs with a thickness of t = 2mm. As a first step to fabricate 
our samples, two-dimensional square arrays of air holes with a lattice constant p — 30mm were drilled through four 
identical copper plates, each hole has the same diameter D = 10mm; Secondly, each air hole was embeded coaxially 
with an annular ring made of teflon with a dielectric constant e r — 2.2. and an appropriate diameter d; Finally, a 
copper cylinder with a diameter of d is embeded in the dielectric ring at the hole center, as shown in the inset of 
Fig.l. Four samples were fabricated with different inner diameters d = 8.3, 6.3, 5.0, 2.0mm respectively. 

The transsmission measurements are performed with Aglient network analyzer 8722ES in a microwave chamber. 
Three pairs of highly directive horn antennas with a gain factor of 24.8dB are utilized to cover the frequency range 
within 5.8^8.2, 8.2~12.4, and!2.4^18.2GHz respectively. The sample plate is vertically positioned on a rotary table, 
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at the center of the line between the emitting and the receiving horns, the horn antennas are 5.0 meters apart from each 
other. Oblique incidence is achieved by rotating the sample plate axially about y axis. Horizontal(i?| \x) or vertical 
polarized (E\\y) Gaussian beams from the horns give rise to the TM or TE polarized incidence. Transmissivity 
throught the sample plate is normalized to the wave energy transmitted between the two horn antennas in free space 
without the sample plate. 



III. MODAL EXPANSION THEROY 

The modal expansion method[l2j is a powerful tool to calculate the transmission of layered system, for example, 
one-dimentional metallic grating. The essence of the method is that the electromagnetic wave fields in the specific 
layers can be decomposed with a series of in-plane eigen-functions correspondingly. The problem can be solved by 
applying boundary continuum condition at all interfaces. At low frequency region, the wave fields inside the aperture 
mainly determined by the fundamental waveguide mode, so that the method can be simplified to one- mode theory [l3T| 
which is fast and convergent. 

With the notions stated above, we develop the modal expansion method to deal with two-dimentional(2D) array 
of MCAAs. We note that the plane wave incidence can not excite TEM guided mode and all the TM guided modes 
in a coaxial aperture accounting for its in-plane symmetry. So we only consider the TE pq guided mode inside the 
MCAAs, where the integer pair (p,q) denotes the mode number [3]. The magnetic field inside the MCAAs within 
the region of < z < h can be written as the superposition of all the forward and backward TE p>q modes: 

&U = E~ j E^ (W^ z - b p , q e^) y), (1) 

where a p g and b p , q are the coefficents of the (p, q) th order of forward and backward guided waves, (3 p>q is the 

wavevector along z axis with (3 p . q — Jk 2 — T^ g ,and T P:9 is the q th root of eigen function J' p {Td/2)N' p (TD/2) — 

J' p {TD/2)N' p (Td/2) — 0. g P , q (x,y) is the in-plane mode profile of the (p,q) th order of guided wave, which can be 
expressed with g p , q (x,y) = gp p>9 (fi,4>)e p + g<f, p>g (p, 4>)h in cylinder coordinates, where 



9 Pp .M <t>) = K ( T P^ d /2) 4 ( T p,iP) ~ 4 ( T P^ d / 2 ) K ( t p*p)] sin W) 

9^M <t>) = - l p IK (Tp-id/ 2 ) J p ( T p-iP) ~ 4 (Tp-id/ 2 ) N p ( T P,,P)] cos(W>) (2) 

The wave fields at the incident side(z < 0) and the outgoing side (z > h)oi the slab can be expanded as the 
superposition of all orders of Bloch waves respectively 

H 1 = I 00 e- in rn e -ik*o, * + y" +0 ° y +0 ° i? mi „ e - l(£ 'i +d '"<"^ie- 4fc ^." z (3) 

' ^ in. — no ' ^ r) — net ' 
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where Iq o represents the coefficent of incident wave, k\\ = kosm(6)e x is the in-plane wavevector, while k ZQ = 
is the wavevector along z axis. R m , n and T m , n are the coefficents of the (m, n)» reflected and transmitted 



Bloch waves. G m ,„ = e x ^p- + Cy^j 1 is the Bloch wavevector of square lattice, where e^and ej,are the unit vectors 
of reciprical space in cartesian coordinate representations. After matching the wave fields of the three regions at two 
interfaces z = Oandz = h, we have 

^ — ^oo ^ — ^oc 

^m,n,p,g(-^0,0^m^n -^m,n) — ^p,q "f" ^p,q 
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k -fTm,n = £°° , E°° , [a P , q e-^ - b p , q e^ h ] (5) 
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where Q m ,n,p,q = fo^dcj) dpe l( - k « +Gm - n ^' r « g^(p, 4>) is an overlap intergral, denoting the coupling coefficent from a 
Bloch wave channel (m, n) th to a guilded wave channel (p, q) th . We expand a plane wave e - lk nP cos (<t>- e ) i n the form of 

Sz^o IT [—iknpcos((j) — 0)] 1 in cylinder coordinates to solve EQs®. Then we derived coefficents of the transmitted 
or reflected Bloch waves for the linear EQs© 

k 

J-m,n — lL m,n,p,q n V u P,9 e u P,q c J 

Pp,q 

k z 

Rm,n — -^0,0^m^n ^m,n,p,q~^ i&p.q ^p,q) (6) 

Pp,q 

As surface resoances are intrinsic properties of the system, we can extract the dispersion of surface resonances by 
assigning a zero value to the incidence, i.e., omitting the term EQQ, and solving the eigen- value 

ofEQs©. 



IV. EVEN-MODE OR ODD-MODE SURFACE RESONANCES 



Circular dots and solid lines in Fig.l) presented the measured and calculated transmission spectra under normal 
incidence (9 = 0°) for samples with difffcrent inner diameters d — 8.3, 6.3, 5.0, 2.0mm. For the sample with d — 8.3mm, 
the measured transmission peak reaches unity at 8.64GHz (Fig.la)). When the inner diameter d becomes smaller, 
the transmission peak at 8.64GHz becomes blue-shifted. In good agreement with the peak and lineshape of the 
measured curves, the calculated transmission spectra also possess two narrow transmission peaks with fano lineshape 
at 9.99GHz and 14.12GHz irrespective to the inner diameter d. These additional transmission peaks are scaled by 
periodicity of the structure. Fig. 2) presents frequencies of the transmission peaks under normal incidence (9 = 0°)with 
respect to different inner diameters d. We calculate the spatial field distributions on the slab to better understand the 
formation and the underlying mechanism of these two kinds of transmission peaks. Fig3) presents the electric field 
component and the power flow component in xz plane on the slab surface with respect to the two transmission peak 
at 8.64GHz and at 9.99GHz. For the narrow resonant peak at 9.99GHz, calculations on the spatial field distribution 
show that the electric field component along z axis is anti-symmetric to the z = ^ plane at both sides of the plate, 
which is rightly the field distribution of an odd-mode surface resonance [Hj], as shown in Fig. 3b) and Fig. 3d). It is 
reasonable that its frequency is mainly decided by the periodicity and symmetry of the array, and indifferent to the 
local geometry of apertures, in that most field energy near the metallic surface distributes far from the apertures with 
the minimum in them, as shown in Fig. 3d). In other words, the odd- mode surface resonance is almost blind to the 
aperture. While the shifted broad transmission peaks in Fig.2) come from the even-mode surface resonant states in 
that the electric field component along z axis is symmetric to the z = § plane at both sides of the plate with most field 
energy localized inside the aperture as well as its vicinity near the surface, as shown in Fig.3a) and Fig.3c) for the 
transmission peak at 8.64GHz respectively. The normalized strength of the electric field inside the aperture can be 
enhanced to 1~2 orders depending on the straitness of the radial gap width (D — d)/2. The field enhancement may be 
applied to the near-field probing or nonlinear optics [l(|. The frequency of an odd-mode surface resonance fm d n{k x ) 

is slightly lower than Rayleigh frequency fm a n lei9h {^/ /) = p y~(Mc + m ) 2 + (k' y + n) 2 where Wood anomaly happens 

[l7l [l8|. Ayy = (k x ,k y ) -p/(27r)is the normalized wavevector. For example, at k x — k y — 0, f^!> lei 9 h ~ Q.QQ3GHz 

is slightly higher than f®Q d — 9.990GHz for the array of holes (d = 8.3mm). With the increase of d, the frequency 

fm d n(k//) suffers an imperceptible blue-shift. All the even- mode branches in Fig.2) are confined within the intervals 

devided by the frequencies {f° d n} °f odd-mode surface resonances, and tends to asymptotic to the lower bound 
with the increase of d. The features imply that the even-mode surface resonances are the collective selections on the 
interplay between the waveguide resonance and in-plane evanescent diffractive channels. 



V. ANGULAR DEPENDENT TRANSMISSION SPECTRA AND DISPERSION DIAGRAM OF 

SURFACE RESONANCES 

We calculate the dispersion diagram of surface resonances for all polarizations, as shown in the solid or dashed 
lines in Fig. 4) for the even-mode or odd-mode surface resonant states.We also measured the transmission spectra 
through all four samples under oblique incidence. The mesured transmissivity is plotted in colormap as a function 
of frequency and the in-plane wavevector fc|| = sm(9)e x , converted from the incident angle 9. For the sample of 
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d = 8.3mm, a transmission peak at / = 8.64GHz in TE polarization is not shifed almost under any incidence angles 
(see FigAa)), reaching unity at normal incidence with 9 = 0°. A similar result has been predicted in a previous study 
on the rectangular holespij. For the other three samples, the first transmission peak in TE polarization become 
dispersive alternatively with respect to the incident angles, the experimental data for d = 6.3, 5.0, 2.0mm is shown in 
Fig Ac, 4e, 4gr) respectively. It implies that the local resonance of the sample with d — 8.3mm lies below the limit 
of Rayleigh frequency while above it for other samples, as we know that a local resonance will couple with different 
Bloch channels when it is below or above Rayleigh frequency. 

Below the limit of Rayleigh frequency j^ yle%g j there only exist G(— 1,0) and G(0,0) Bloch channels. Without 
the projection of wave vector component of the incidence along y axis, the TE polarized waves incident in the y = 
plane can not couple with these Bloch channels. So for d — 8.3mm, without the assistance from any Bloch channels, 
we can observe a flat polariton branch in TE mode, as shown in Fig. 4a). Meanwhile the TM polarized incidence in 
the y = plane can be coupled to the G(— 1,0) Bloch channel due to the same reason; the TM poliarized surface 
resonant states in even mode comes from the interference between the waveguide resonance and the G(— 1,0) Bloch 
channels, giving rise to two anti-symmetric asymptotic branches in FigAb). The bright curvatures for measured data 
prove the existance of the two branches, as shown FigAb). Thus the role of local resonance on the surface waves 
below JiQ is revealed. 

More diffractive channels lie above f^ yei3 , such as G(0,±1), andG(±l,±l) etc. When d becomes smaller, the 
waveguide resonance will shift across f^y lei 9 h f rom below. TE polarized incidence will couple with the waveguide 
resonance via G(0, ±1), andG(— 1, ±1) channels, resulting in different features of dispersion diagrams and transmission 
spectra. As an example, for the sample with d = 6.3mm, the three branches below fif d = 14.12GHz are rightly the 

consequence of the interplay between the local resonance and the G(0, ±1), G(— 1, ±1) channels under TE polarized 
incidence. As shown in Fig. 4c), the bright-colored curvature exactly fits the curve of the lowest branch with the 
brightest point at the location (k x = ir/p, f = 9.97GHz), which is the four-fold degenerate point about G(— 1, ±1), 
and G(0, ±1). The maximal transmissivity at the Brillouin zone boundary means that the strongest coupling between 
the local resonance and the Bloch channels is likely to occur at the degenrate states with high symmetry. The 
coupling between local resonance and diffractive channels will induce highly inhomogeneous field distribution on the 
surface, resulting in the opening of bandgaps, similar to the case of corrugated metallic surface 19]. As shown in 
Fig. 4c), existence of two bandgaps between the three branches verifies the analysis. The coupling will inevitably 
induce the split of the frequency at the degenrate points. We notice that right above the "brightest" state (k x — n/p, 
f = 9.97GHz) , the second TE polarized branch in even mode converge with four odd-mode branches (—1, ±1) and 
(0, ±1) at the degenerate state(fc x = n/p, f =11.12GHz) on BZ boundary. The degeneracy demonstrates that the 
state also possesses anti-symmetric distribution of the spatial field and minimum transmissivity. It also applies for 
Fig.4e) and Fig.4g); the superposed bright curvature on the second branch becomes weaker, reaching minimum at 
the four-fold degenerate state. The rule governs all dispersion diagrams in TM polarization, on the second branch 
at k x = ir/p, the third branch at k x = n/2p and k x = 3ir/2p, the fourth branch at k x = ixjlp in Fig.4b,d,f,h). 
Transmission extrema at the degenerate states are understandable. Considering these brightest states with four-fold 
degeneracy, they are likely to couple with the incident waves with strongly localized field within the apertures and 
have zero group velocity along any direction within the xy plane. Meanwhile, the odd-mode surface resonant states 
are invisible in measurements, and the geometry of apertures has little impact on the odd-mode branches. 



VI. CONCLUSION 



In summary, we perform microwave experiments and rigorous modal expansion method to investigate the enhanced 
microwave transmission through the array of metallic coaxial annular apertures (MCAAs) as well as the dispersion 
of surface resonances. The even-mode and the odd-mode surface resonances are clarified theoretically by the spatial 
field distributions and the dispersion diagram. The impact of local resonance is thoroughly embodied in the even- 
mode surface resonant states. While the odd-mode surface resonant states, invisible in measurements, are scaled by 
periodicity, invariant to different local geometry of the unit cell. The enhanced transmission is the collective selections 
on the interplay between the local resonances and the evanescent Bloch channels on the surface. Transmissivity 
reaches extrema at the degenerate points in the dispersion diagram of surface resonances. The work provides a 
holisitc description on the interplay between the local resonance and the plane wave excitations via evanesent Bloch 
channels. 
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Figure captions: 

FIG.l.Simulated(lines) and Measured(open circles) transmission spectra through the array of MCAAs with 
different inner diameter d = 8.3mm(a),6.3mm(b), 5.0mm(c), 2.0mm(d) respectively. The inset presents the photo of 
the sample (d = 8.3mm) 

FIG. 2. Resonant frequencies of even mode(solid lines) and odd mode(dashed lines) under normal incidence (fcii = 0) 
with respect to the inner diameter d varying in the range of 0.0 ~ 4.9mm. The measured frequencies for different 
samples are marked with stars. 

FIG. 3. Spatial field distributions of the even/odd surface mode in xz plane for the sample with d — 8.3mm. (a) E 
field, even; (b) E field, odd; (c) Power flow, even; (d)Power flow, odd 

FIG. 4. Dispersion relations of even surface resonances (solid lines) and odd surface surface resonances (dashed lines) 
for samples with different inner diameter d = 8.3, 6.3, 5.0, and 2.0mm. (a), ( c), (e), (g) for TE polarization and (b), 
(d),( f), (h) for TM polarization. Frequency / is normalized to the limit of Rayleigh frequency cO/p = 9.993GHz. 
Measured transmittivities with respective to the in-plane wave vector and frequency are plotted in color-scale. 
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FIG.l.Simulatcd(lines) and Measured(open circles) transmission spectra through the array of MCAAs with different 
inner diameter d = 8.3mm(a),6.3mm(b), 5.0mm(c), 2.0mm(d) respectively. The inset presents the photo of the 

sample (d = 8.3mm) 
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FIG. 2. Resonant frequencies of even mode(solid lines) and odd mode(dashed lines) under normal incidence (fcii = 0) 
with respect to the inner diameter d varying in the range of 0.0 ~ 4.9mm. The measured frequencies for different 

samples are marked with stars. 
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FIG. 3. Spatial field distributions of the even/odd surface mode in xz plane for the sample with d = 8.3mm. (a) E 
field, even; (b) E field, odd; (c) Power flow, even; (d)Power flow, odd 
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FIG. 4. Dispersion relations of even surface resonances (solid lines) and odd surface surface resonances (dashed lines) 
for samples with different inner diameter d = 8.3, 6.3, 5.0, and 2.0mm. (a), ( c), (e), (g) for TE polarization and (b), 
(d),( f), (h) for TM polarization. Frequency / is normalized to the limit of Rayleigh frequency cO/p — 9.993GHz. 
Measured transmittivities with respective to the in-plane wave vector and frequency are plotted in color-scale. 



